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ABSTRACT 

We present high-quality Keck/LRIS longslit spectroscopy of a pilot sample of 25 local active galaxies 
selected from the SDSS (0.02 < z < 0.1; Mbh> lO^M©) to study the relations between black hole 
mass (Mbh) and host-galaxy properties. We determine stellar kinematics of the host galaxy, deriving 
stellar- velocity dispersion profiles and rotation curves from three spectral regions (including CaH&K, 
Mgib triplet, and Call triplet). In addition, we perform surface photometry on SDSS images, using a 
newly developed code for joint multi-band analysis. BH masses are estimated from the width of the H/3 
emission line and the host-galaxy free 5100A AGN luminosity. Combining results from spectroscopy 
and imaging allows us to study four Mbh scaling relations: Mbh-ct, MBH-isph, MBH-Mgph,*, and 
MBH-M3ph,dyn. We find the following results. First, stellar-velocity dispersions determined from 
aperture spectra (e.g. SDSS fiber spectra or unresolved data from distant galaxies) can be biased, 
depending on aperture size, AGN contamination, and host-galaxy morphology. However, such a bias 
cannot explain the offset seen in the Mbh-c relation at higher redshifts. Second, while the CaT 
region is the cleanest to determine stellar-velocity dispersions, both the MgIb region, corrected for 
Fell emission, and the CaHK region, although often swamped by the AGN powerlaw continuum and 
emission lines, can give results accurate to within a few percent. Third, the Mbh scaling relations of 
our pilot sample agree in slope and scatter with those of other local active and inactive galaxies. In 
the next papers of the series we will quantify the scaling relations, exploiting the full sample of ^ 100 
objects. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: active — galaxies: 
evolution — quasars: general 



1. INTRODUCTION 

The empirical relations between the mass of 
the central supermassive black hole (BH) and the 
properties of the spheroid (ellipticals and classical 
bulges of spirals) such as stellar-velocity disper sion a 
(jFerrarese fc Merrittl 120001: iGebhardt et al.l[2000[ l. stel- 
lar mass fe.g.. iMarconi fc HuntI 120031). and luminos- 
ity (e.g., iHaring fc Rixl 1200 ^ ^discovered in the lo- 
cal Universe have been interpreted as an indica- 
tion of a close connection between the growth of 
the BH and the formation and evo l ution of galax- 

"20001: iVolonteri et al.l 



ies (e . g., |K ai 
200a ICiotti 



K auffma nn fc Haehneltl 
fc Ostrikei * '2007!: 



Di Matteo et al.l 



, , , _ _ .Hopkins et al.l l2007t 

20081 : iHopkins et al.l |2009() . In this 



framework, active galactic nuclei (AGNs) are thought to 
represent a stage in the evolution of galaxies in which the 
supermassive BH is actively growing through accretion. 

To understand the origin of the BH mass (Mbh) scal- 
ing relations, our gro up has been stud ying their e volu- 
tion with cosmic time dJ reu et al."2004':' Woo et al.ir2006l: 
iTreu et al.l [20071: iWoo e t al. 2008; Bcnne rt et al.l I2010D . 
To distinguish mechanisms causing evolution in a (e.g., 
dissipational merger events) and L^-ph (e.g., through pas- 
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sive evolution due to aging of the stellar population, or 
dissipationless mergers), we simultaneously study both 
the Mbh-c and MBH-^sph relations for a sample of low- 
luminous AGNs, Seyfert-1 galaxies, at z ~ 0.4 and 0.6 
(lookback time 4-6 Gyrs). Our results reveal an off- 
set with respect to the local relationships which can- 
not be accounted for by known systematic uncertain- 
ties. The evolutionary trend we find (e.g., Mbh/^spIiOC 
jjj_^^-ji.4±o.2^ including selection effects; iBennert et al.l 
l2010f ) suggests that BH growth precedes spheroid assem- 
bly. Several other studies have found results in qualita- 
tive agreement with ours, over different ranges in black 
hole mass and redshifts, and with different observing 
techniques (eg.. I Walter et all [20 04: Shields ct al. 200a 
McLure et all l2006t iPeng et al.lT2006bi: ,Salviander et all 



20071: iWeiss et al.r [2007t iRiechers et all 120081 



2009; 



Gu et al. ''2009': 'Ja hnke et al.1 120101 : iDecarh et'SlboiC : 
, Merloni et al. 2010). 

However, to study the evolution of the Mbh scaling 
relations, it is crucial to understand slope and scatter 
of the local relations. In particular, an open question 
is whether quiescent galaxies and active galaxies follow 
the same relations, as expected if the nuclear activity 
was just a tra nsient phas e in t he life-cycle of galax- 
ies. Recently, iWoo et al] (|2010[) presented the Mbh- 
cr relation for a sample of 24 active galaxies in the lo- 
cal Universe, for which the BH m ass was derived via 
reverberation mapping : (RM) (e.g., IWandel et all [1991 



iKaspi et al.|[2000ll2005l: iBentz et al.ll2006l l2009bi r They 
find a slope (/3 = 3.55 ± 0.60) and intrinsic scatter 
(cint = 0.43 ± 0.08) which are indistinguishable from 
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that of quiescent ga laxies fe.g.. iFerrarese fc FordI [20051 : 
iGiiltekin et al.l [20091) within the uncertainties, support- 
ing the scenario in which active galaxies are an evolu- 
tionary stage in the life cycle of galaxies. 

While the great advantage of such a study is the multi- 
epoch data which provide more reliable measurements of 
the BH mass, such a quality comes at the expense of 
quantity. Studies based on larger samples drawn from 
the Sloan Digital Sky Survey (SDSS) infer the BH mass 
indirec tly from single-ep och spectra fe.g.. iGreene fc Hoi 
I2006at iShen et al.ll2008[ ). They hinted at a shallower 
AfeH-c relation than that observed for quiescent sam- 
ples, but the available dynamic range is too small to be 
conclusive. In particular, the relation above 10^'^ Mq is 
very poorly known, which has profound implications for 
evolutionary studies that by necessity focus on this mass 
range. 

However, there is another uncertainty in the Mbh- 
(T relation that arises when measuring a from fiber- 
based SDSS data (e.g., IGreene fc H"oll2006at iShen et al.l 
120081 ) and also from the unresolved "aperture spectra" 
for more distant galaxies (e.g., our studies on the evolu- 
tion of the Mbh-ct relation; Wo o et a l. 2006, 2008, J.-H. 
Woo et al. 2011, in preparation). Local active galaxies 
seeni to span a range of morphologies (e.g.,lMalka n et ah 
1998tlHunt fc Malkan|[200llKim et al . 2008; Bcnt z et al. 
2009al) and a significant fraction (> 15/40) of our dis- 
tant sample of Seyfert-1 galaxi es have morphologies of 
Sa or later (jBennert et al.ll2010t) . Given the diversity of 
morphologies of AGN hosts, it is most likely that there 
is a degree of rotational support: If the disk is seen 
edge-on, the disk rotation can bias a towards higher val- 
ues. However, since the disk is kinematically cold, it can 
also result in the opposite effect, i.e. biassing a towards 
smaller values, if the disk is seen face on (e.g.. Woo ct al] 
I2006f ). Either way, it questions the connection between 
the "global" dispersion measured by those experiments 
and the spheroid-only dispersion which may in fact scale 
more tightly with BH mass. 

More generally, measuring a in type-1 AGNs is com- 
plicated by the presence of strong emission lines and a 
continuum that dilutes the starlight, a can be measured 
from different spectral regions with different merits and 
complic ations ([Greene fc Holl2006bl for inactive galaxies 
see also iBarth et al.ll2002f ). Finally, there is are differ- 
ent (T measurements in use: e.g., the luminosity- weighted 
line-of-sight velocity dispe rsion within the sphe r oid ef - 
fective radius (crcfi, e.g., IGebhardt etHI [200(1 [200l . 
and the central velocity dispersions normalized to an 
aperture of radius equal t o 1/8 of the galaxy ef- 
fective radius (cri/srcff: e.g.. IFerrarese fc Merritti [20001 : 
IFerrarese fclbr d 2005). 

Shedding light on the issues outlined above is essen- 
tial to understand what aspects of galaxy formation and 
AGN activity are connected, but it requires spatially- 
resolved kinematic information for a large sample of local 
AGNs. For this purpose, we selected a sample of ^100 
local (0.02 < 2 < 0.1) Seyfert-1 galaxies from the SDSS 
(DR6) with Mbh> lO^M© and obtained high-quality 
longslit spectra with Keck/LRIS. From the Keck spec- 
tra, we derive the BH mass and measure the spatially- 
resolved stellar-velocity dispersion from three differ- 
ent spectral regions (around CaH+K AA3969,3934A; 



around Mgib AAA5167,5173,5184A triplet; and around 
Ca IIAAA8498,8542,8662A triplet). The spectra are com- 
plemented by archival SDSS images (g', r', i', z') on which 
we perform surface photometry using a newly developed 
code to determine the spheroid effective radius, spheroid 
luminosity, and the host-galaxy free 5100A luminosity 
of the AGN (for an accurate BH mass measurement). 
Our code allows a joint multi-band analysis to disentan- 
gle the AGN which dominates in the blue from the host 
galaxy that dominates in the red. The resulting multi- 
filter spheroid luminosities allow us to estimate spheroid 
stellar masses. 

Combining the results from spectroscopic and imaging 
analysis, we can study four different BH mass scaling 
relations (namely Mbh-ct, A/eH-isph, MBK-Msph,*, and 
AfBH--^^sph,dyn)- In this paper, we focus on the methodol- 
ogy and present results for a pilot sample of 25 objects. 
The full sample will be discussed in the upcoming pa- 
pers of this series. The paper is organized as follows. 
We summarize sample selection and sample properties 
in § [21 observations and data reduction in § [3l § [H de- 
scribes the derived quantities, such as spatially-resolved 
stellar-velocity dispersion and velocity, aperture stellar- 
velocity dispersion, BH masses, surface photometry, and 
spheroid masses. In §[5l we describe comparison samples 
drawn from literature, consisting of local inactive and ac- 
tive galaxies. We present and discuss our results in § [6l 
including the BH mass scaling relations. We conclude 
with a summary in § [7] In Appendix [3 we describe the 
details of a python-based code developed by us to deter- 
mine surface-photometry from multi-filter SDSS images. 
Throughout the paper, we assume a Hubble constant of 
Ho = 70kms-i Mpc-\ Qa = 0.7 and = 0.3. 

2. SAMPLE SELECTION 

Making use of the SDSS DR6 data release, we selected 
type-1 AGNs with Mbh> 1O''M0, as estimated from the 
spectra based on their optical lum inosity and H/3 fuU- 
width-at-half-maximum (FWHM) (jMcGih et al.l 120081 ). 
We restricted the redshift range to 0.02 — 0.1 to ensure 
that both the Ca triplet and a bluer wavelength region 
are accessible to measure stellar kinematics and that the 
objects are well resolved. This results in a list of 332 ob- 
jects from which targets were selected based on visibility 
during the assigned Keck observing time. Moreover, we 
visually inspected all spectra to make sure that the BH 
mass measurement is reliable and that there are no spu- 
rious outliers lacking broad emission lines (~5% of the 
objects). A total of 111 objects were observed with Keck 
between January 2009 and March 2010. Here, we present 
the methodology of our approach and the results for our 
pilot sample of 25 objects. Their properties are summa- 
rized in Table [TJ Fig. [T] shows postage stamp SDSS-DR7 
images. The results for the full sample will be presented 
in the forthcoming papers of this series. 

All 25 objects were covered by the VLA FIRST (Faint 
Images of the Radio Sky at Twenty-cm) survey^, but only 
10 have detected counterparts within a r adius of 5". Out 
of these 10, seven are listed in lLi et all (|2008[ ) with only 
one being radio-loud. Thus, the majority (>~85%) of 
our objects are radio-quiet. Note that none of the ob- 

^ See VizieR Online Data Catalog, 8071 (R. H. Becker et al., 
2003) 
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jects has HST images available. As our sample was se- 
lected from the SDSS, most objects are inc luded in stud- 
ies th at discuss the local BH mas s function (iGreene 
I2007| ) or BH fundamental plane (|Li et al.ll2008D In addi 
tion, 1535+5754 (Mrk 290) has a reve rberation-mapped 
BH mass from iDenney et all (|2010t ). We wih com- 
pare the BH masses derived in these studies with ours 
when we present the full sample. For a total of eight 
objects, stellar-velocity dispersion measurements from 
aperture spectra exist in the l iterature (m ainly derived 
from SDSS fiber data: six in IGreene fc Ho 2006a and 
five in lShen et al. 2008', with thre e overlapping, and one 
object in lNeison fc Whittlelll996l determined from inde- 
pendent spectra, but included in both SDSS studies). We 
briefly compare the stellar- velocity dispersions derived in 
these studies with ours in § 14.21 but will get back to it in 
more detail when we present the full sample. 

3. SPECTROSCOPY: OBSERVATIONS AND DATA 
REDUCTION 

We here summarize only the spectroscopic observations 
and data reduction. The photometric data consist of 
SDSS archival images and the data reduction is summa- 
rized in Appendix lAl 

All objects were observed with the Low Resolution 
Imaging Spectrometer (LRIS) at Keck I, using a 1" wide 
longslit, the D560 dichroic, the 600/4000 grism in the 
blue, and the 831/8200 grating in the red (central 
wavelength = 8950A). In addition to inferring the BH 
mass from the width of the broad H/3 line, this setup 
allows us to simultaneously cover three spectral re- 
gions commonly used to determine stellar-velocity dis- 
persions. In the blue, we cover the region around the 
CaR+K AA3969,3934A (hereafter CaHK) and around the 
Mgib AAA5167,5173,5184A triplet (hereafter Mgib); in 
the red, we cover the Ca IIAAA8498,8542,8662A triplet 
(hereafter CaT). The instrumental spectral resolution is 
~90 kms~^ in the blue and '~45 kms~^ in the red. 

The long slit was aligned with the host galaxy major 
axis as determined from SDSS ("expPhi_r"), allowing us 
to measure the stellar-velocity dispersion profile and ro- 
tation curves. Observations were carried out on January 
21 2009 (clear, seeing 1-1.5"), January 22 2009 (clear, 
seeing ~1.1"), April 15 2009 (scattered clouds, seeing 
'-I"), and April 16 2009 (clear, seeing -0.8"; see also 
Table 

Note that all data included in this paper were obtained 
before the LRIS red upgrade. The rest of our sample 
(—75 objects) benefited from the upgrade with higher 
throughput and lower fringing (data obtained from June 
2009 onwards) and will be presented in an upcoming pa- 
per (C. E. Harris et al. 2011, in preparation). A total of 
25 objects were observed with the old red LRIS chip. For 
three objects, the spectra did not allow a robust measure- 
ment of the stellar kinematics, due to dominating AGN 
continuum and emission lines. However, we were able to 
determine BH mass and surface photometry. Thus our 
"imaging" sample consists of 25 objects, our "spectro- 
scopic" sample of 22 objects. 

The data were reduced using a python-based script 
which includes the standard reduction steps such as bias 
subtraction, fiat fielding, and cosmic ray rejection. Arc- 
lamps were used for wavelength calibration in the blue 



spectral range, sky emission lines in the red. AOV Hip- 
parcos stars, observed immediately after a group of ob- 
jects close in coordinates (to minimize overhead), were 
used to correct for telluric absorption and perform rela- 
tive flux calibration. 

From these final reduced 2D spectra, we extracted ID 
spectra in the following manner. For the blue, a central 
spectrum with a width of 1.08" (8 pixels) was extracted 
to measure the H/3 width for BH mass determination, 
i.e. encompassing the broad-line region (BLR) emission 
given a typical seeing of ~1" and a slit width of 1". To 
measure the stellar-velocity dispersion and its variation 
as a function of radius, we extracted a central spectrum 
with a width of 0.54" (0.43") for the blue (red). Outer 
spectra were extracted by stepping out in both direc- 
tions, increasing the width of the extraction window by 
one pixel at each step (above and below the trace) choos- 
ing the stepsize such that there is no overlap with the 
previous spectrum. If the S/N of the extracted spectrum 
fell below 10 pix~^ (at rest wavelength 5050-5450A in 
the blue, and 8480-8690A in the red), the width of the 
extraction window was increased until an S/N of at least 
10 pix~-^ was reached. We only use spectra with S/N 
> 10 pix^^. The signal-to-noise ratio (S/N) of the final 
reduced total spectra (extraction with aperture radius 
of —7") is on average —80 pix~^ in the blue (ranging 
from —30 pix^^ to 190 pix~^) and —70 pix~^ in the red 
(ranging from —20 pix~^ to 170 pix~^; Table [2]). 

The majority of objects (15/22) display broad nuclear 
Fell emission in their spectra (— 5150-5350A), compli- 
cating the fj measurements from the Mgib region. For 
those objects, we fitted a set of IZwl templates, with 
various widths and strengths, in addition to a feature- 
less AGN continuum. The best fit was determined by 
minimizing and then subtracted. Details of this pro- 
cedure are given in iWoo et all pool . We first derived 
the best fit for the central spectrum and then used the 
same Fell width also for the two outer spectra that 
are still affected by Fell due to seeing effects and slit 
width. In Fig. [21 we show an example of the Fell sub- 
traction. Fig. [3] compares the observed spectra to the 
Fell emission-subtracted spectra. Note that we did not 
correct the CaHK region of our spectrum for Fell emis- 
sion, since th e broad Fell features near —3950 are weaker 
and broader (I Greene infoll2006Hl . 

Details of the spectroscopic observations and data re- 
duction are given in Table [2l 

4. DERIVED QUANTITIES 

In this section, we summarize the results we derive 
both from the spectral analysis (stellar-velocity disper- 
sion, velocity, and H/3 width) and image analysis (surface 
photometry and stellar masses) as well as from combin- 
ing results from both (BH mass and dynamical masses). 
We will also distinguish between different stellar-velocity 
dispersion measurements and define the nomenclature we 
use. 

4.1. Spatially- Resolved Stellar- Velocity Dispersion And 

Velocity 

The extracted spatially-resolved Keck spectra allow us 
to determine the stellar-velocity dispersion as a func- 
tion of distance from the center. In the following, we 
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TABLE 1 

Sample Properties 
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Note. — Col. (1): Target ID used throughout the text (based on RA and DEC). Col. (2): Full SDSS name. Col. (3): Rcdshift from SDSS-DR7. 
Col. (4): Luminosity distance in Mpc, based on redshift and the adopted cosmology. Col. (5): Scale in kpc/arcsec, based on redshift and the 
adopted cosmology. Col. (6): Right Ascension. Col. (7): Declination. Col. (8): i' AB magnitude from SDSS-DR7 photometry ( "modelMagJ" ) . 
Col. (9): Alternative name(s) from the NASA/IPAC Extragalactic Database (NED). 



will refer to these measurements as <Tspat) in contrast 
to velocity dispersion determined from aperture spec- 
tra, as discussed in the next subsection. The advan- 
tage of spatially-resolved spectra is twofold: For one, 
the spatially-resolved stellar velocity dispersions are not 
broadened by a rotating disk (if seen edge-on) and sec- 
ond, the contamination by the AGN powerlaw contin- 
uum and broad emission lines will only affect the nuclear 
spectra, but not spectra extracted further out. 

A python-based code measures the stellar-velocity dis- 
persion from the extracted spectra, using a linear com- 
bination of G&K giant stars taken from the Indo-US 
survey, broadened to a width ranging between 30-500 
kms"'^. In addition, a polynomial continuum of order 3- 
5 was fitted, depending on the object and fitting region. 
The code uses a Markov-chain Monte Carlo (MCMC) 
simulatio n to find th e best -fit velocity dispersion and 
error; see ISuvu et all (pOlOl) for more detailed descrip- 
tion of the fitting technique. Three different regions 
were fitted: (i) the region around CaT, 8480-8690A; 
(ii) the region around Mgib that also includes several 
Fe absorption lines, 5050-5450A (i.e. redwards of the 
[OIII] A5006.85A and the HeIAA5015.8,5047.47A AGN 
emission lines); and (iii) the region around the CaH-|-K 
lines, 3735-4300A (i.e. bluewards of the broad H7 and 
[OIII] A4363.2IA AGN emission lines). In Fig. HlTl ex- 
amples of fits are shown for all three regions. 

In region (i), the broad AGN emission line 
OIA8446A contaminates the spectra in the central re- 
gion for some objects. In those cases, we excluded the 
first CaT line and only used the region 8520-8690. In 



objects at higher redshifts, the third CaT line can be af- 
fected by telluric absorption (although we attempted to 
correct for this effect), and had to be excluded in some 
cases. 

In region (ii), several AGN narrow-emission lines 
were masked, if present, such as [Fe VI] AAA5I45.77, 
5176.43, 5335.23A; [Fe VII] AA5158.98, 5277.67A; 
[NI] AA5197.94,520 0.4lA: and f CaV] A5309.1 8A (wave- 
lengths taken from lMoor(j|1945l: IBowenlll960D . The blue 
spectra end at an observed wavelength of ~5600A, which 
corresponds to restframe ~5200A for our highest red- 
shifted object (z=0.076). 

Region (iii) includes AGN emission from e.g., [FeV] 
A4227.49A and various broad Hel and Balmer lines, that 
were masked during the fitting procedure. Also, the 
CaHA3969A hue is often filled by AGN emission (He). 
Thus, in most cases, we restricted the fitting region to 
4150-4300A for the central spectra (see Fig. [B]), due to 
AGN contamination. Only in the outer parts, the wave- 
length regime 3735-4300A was fitted (Fig. [T]). In the 
following, we still generally refer to region (iii) as CaHK 
region, even though it might not actually include the 
CaH+K line in the central spectra where it is contami- 
nated by AGN emission. 

The code used to determine the stellar- velocity disper- 
sion also gives the line-of-sight velocity. (Note that we 
set the central velocity to zero.) The resulting measure- 
ments for both stellar-velocity dispersion as well as ve- 
locity as a function of distance from the center are shown 
in Figures mUni 

The error bars are often higher in the center due to 
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Fig. 1. — Postage stamp SDSS-DR7 multi-filter images. North is up, east is to the left. The size of the field-of-view is 50"x50" (cor- 
responding to ~20x20kpc up to ~76x76kpc for the redshift range covered by the objects presented here). For 0206-0017, we show the 
central 100" X 100". 

AGN contamination. This is also the reason why the 
error bars for velocity dispersions determined from the 
Mgib or CaHK region arc higher, since in the blue, the 
contamination by the AGN powerlaw continuum and 
broad emission lines is more severe than for the CaT 
region. On the other hand, in the outermost spectra, the 
S/N is the dominating error source. 

For comparison with literature, we choose the ve- 
locity dispersion determined from the CaT region as 
our benchmark, since th is is the region le ast affected 
by template mismatch (jBarth et al.l 1200^ and AGN 
contamination fro m a featureless con tinuum as well 
as emission lines (IG reene "fcHol[2006B . We calculate 
the luminosity-weighted line-of-sight velocity dispersion 
within the spheroid effective radius (determined from the 
surface photometry of the SDSS DR7 images as outlined 
below): 



Fig. 2. — Example of the broad Fell emission fit using a IZwl 
template and a featureless continuum for 0121-0102 (central spec- 
trum; ~0.54"xl"). This object has the strongest Fell emission 
in our sample. Upper panel: Observed spectra with best-fit 
model (dotted line) composed of a broadened Fell template and 
AGN+stellar continuum (dashed line). Lower panel: Broadened 
Fell template. 
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TABLE 2 

Details of Observation and Reduction 



Object 


PA 


Exp. time 
(sec) 


Date 


S/Nblue 


S/N,,d 


Fell sub. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


ni oi m no 


65.6 


1200 


m oi nn 


111.7 


85.8 


yes 


UiUb-UUi ( 


176.0 


1200 


ni oo nn 


152.7 


142.8 


no 




171.2 


1200 


m oo nn 

ui-z^-(jy 


50.2 


40.2 


yes 


U8U2+alU4 


82.9 


1200 


ni oi nn 

ui-2i-uy 


79.1 


72.8 


yes 


Uo4d+zozz 


50.9 


1200 


ni oo nn 

ui-zz-uy 


103.6 


89.2 


no 


lU4z+U4i4 


126.2 


1200 


n /I T a nn 
U4-lD-Uy 


53.8 


52.1 


yes 


lU4ci+ilUo 


128.2 


600 


r\ A An nn 


22.9 


17.4 


no 


iU4y+z40i 


29.9 


600 


n /I 1 nn 
U4-iD-Uy 


52.2 


y1 Q 1 
43.1 


yes 


1 1 m 1 1 1 no 


147.5 


600 


n /I T nn 
U4-lD-Uy 


32.2 


37.5 


yes 


T 1 1 C 1 /1100 

lilD+41z.3 


11.7 


850 


c\ A "1 CT nn 

04-15-09 


46.3 


62.4 


no 


1144+0000 


20.7 


600 


04-16-09 


58.7 


62.4 


no 


IzlU+oozU 


0.8 


600 


04-16-09 


113.1 


132.7 


yes 


lZ0U-Uz4y 


73.9 


1200 


n /I 1 nn 
04-16-09 


37.3 


45.8 


yes 


1323+2701 


8.1 


700 


04-16-09 


25.0 


28.7 


no 


1355+3834" 


78.0 


300 


04-16-09 


34.7 


34.4 


no 


1405-0259 


64.8 


1600 


04-16-09 


54.4 


72.2 


yes 


1419+0754 


19.3 


900 


04-16-09 


58.8 


80.1 


yes 


1434+4839 


152.1 


600 


04-16-09 


49.6 


61.3 


yes 


1535+5754 


103.8 


1200 


04-15-09 


180.8 


169.6 


yes 


1545+1709 


60.0 


1200 


04-15-09 


83.7 


91.3 


no 


1554+3238 


169.1 


1200 


04-15-09 


83.6 


95.0 


yes 


1557+0830" 


58.6 


1200 


04-15-09 


54.9 


55.6 


yes 


1605+3305 


90.2 


1200 


04-15-09 


80.0 


80.9 


yes 


1606+3324 


20.8 


1200 


04-15-09 


44.9 


59.1 


yes 


1611+5211 


114.3 


1200 


04-15-09 


72.3 


84.4 


no 



Note. — Col. (1): Target ID. Col. (2): Position angle of major axis, along 
which the long slit was placed (taken from SDSS-DR7 "expPhi_r"). Col. (3): Total 
exposure time in seconds. Col. (4): Date of Observations ( month-day- year) . Col. 
(5): S/N in total blue spectrum (per pix, at rest wavelength 5050-5450A), aperture 
size ~7". Col. (6): S/N in total red spectrum (per pix, at rest wavelength 8480- 
8690A), aperture size ~7". Col. (7): Subtraction of broad Fell emission. 

Note that for these three objects, the spectra did not allow a robust measurement 
of the stellar kinematics due to AGN contamination and we only present BH mass 
and results from surface photometry. 

with Ur) = I jreS) ■ exp(-7.67 • [(r/rrefi)^-^^ - 1]) 
Ide VaucouleuTsI (pisl ) profile. (The range "-reff" to 
"-|-reff" refers to the fact that we extracted spectra sym- 
metrically around the center of each object, along the 
major axis, and measured stellar velocity dispersions 
from each of them; see also Figs.lSlfTUl) As our CTgpat mea- 
surements are discrete, we interpolate over the appropri- 
ate radial range using a spline-function. In the follow- 
ing, we refer to the spatially-resolved stellar-velocity dis- 
persion within the spheroid effective radius as crspat,rcff- 
This represents the spheroid-only dispersion within the 
effective radius, free from broadening due to a rotating 
disk component. Note that the only place where we 
show spatially-resolved velocity dispersion at a certain 
radius is in Figures [SlfTUl (crgpat) and in Figure [T5b for 
the ratios; otherwise we always refer to the luminosity- 
weighted spatially-resolved velocity dispersion within the 
effective radius Uspat.rcfr- 

We can in principle choose arbitrary integration limits, 
to e.g. determine the stellar-velocity dispersion within 
one-eighth of the effective radius, one-half of the effec- 
tive radius, or 1.5" (comparable to SDSS fiber spectra). 
While these are all "radii" found in the literature, we 
will not be using them for comparison in this paper, 
since literature values all refer to aperture data (and not 
spatially-resolved as discussed here; see next section). 

4.2. Aperture Stellar- Velocity Dispersion 



While we consider CTgpat.reff the spheroid dispersion 
within the spheroid effective radius free of disk contam- 
ination, comparison with literature data (such as fiber- 
based SDSS data or unresolved aperture spectra for more 
distant galaxies) requires us to also determine aperture 
stellar-velocity dispersions. To do so, we extracted ID 
spectra by increasing the width of the extraction win- 
dow by one pixel, leaving the centroid fixed to the central 
pixel. We then measured the stellar-velocity dispersion 
for each of the extracted spectra using the same proce- 
dure as outlined above. We use the same mask/fitting 
region as used for the center for the spatially-resolved 
spectra (since here, all spectra will suffer from the AGN 
contamination). Also note that broad Fell emission was 
subtracted for all rows, if present, fixed to the widths de- 
termined in the center. We refer to the resulting stellar- 
velocity dispersion from these spectra as "aperture" CTap. 
Note that Cap not only contains the AGN contamination 
(continuum and emission lines) independent of extracted 
aperture (and that is why here the fitting region that 
we refer to as CaHK region does not actually include the 
CaH+K line; see § l4.ip but can also be broadened by any 
rotational component present or biased to lower values 
in case of contribution of a kinematically cold disk seen 
face on. Also, the resulting fTap,rcff is already luminosity 
weighted due to the way the spectra are extracted. The 
results are shown in Fig. 1111 
We determine an aperture stellar-velocity dispersion 
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Fig. 3. — Central spectra (~0.54"xl"; normalized to average flux over plotted wavelength region), before (gray) and after (black; shifted 
arbitrarily by -0.2 for comparison) broad Fell subtraction for the 15 objects for which Fell was subtracted. The location of the stellar 
absorption lines MgIbA5175A and FeA5270A is indicated by dotted lines. While some objects show prominent stellar absorption features, 
for others the region is swamped by AGN emission. Note that the spectra end at the rest wavelength 5600A/(l+z). 



within the effective radius (7ap,rGff by choosing the aper- 
ture size identical to the spheroid effective radius of a 
given object. 

To compare our resuhs with SDSS fiber spectra, we de- 
termine (Tap, 1.5", measured from aperture spectra within 
the central 1.5" radius as a proxy for what would have 
been measured with the 3" diameter Sloan fiber. Note, 
however, that in fact, our (Jap, 1.5" corresponds to a rect- 
angular region with 1.5" radius and 1" width, given the 
width of the long slit. For eight objects, we can di- 
rectly compare our results for Jan. 1 .."s" with thos e de- 
rived from SPSS fibe r spectra bv lShen et all (|2008[ ) and 
I Greene fc Hoi (j2006a[ ). While individual objects can dif- 
fer by up to 25%, slightly larger than the quoted uncer- 
tainties for the SDSS spectra (~10-15%), on average, the 
measurements agree very well ((Tap 1 5"/(tsdss = 1-05 ± 

We summarize the different stellar- velocity dispersion 



measurements in Table |3l for the CaT region only.^ 



4.3. Black Hole Mass 

Black hole masses are estimated using the 
empirically calibrated photo-ionization method, 
also sometimes kn o wn as "v irial me thod" (e.g., 
Wandel et ai] [1991 iVe stergaard 2002*: iWoo fc Urrvl 



200a IVestergaard fc Peterson. 2006.: .McGill et al.ll2008[) . 



Briefly, the method assumes that the kinematics of the 
gaseous region in the immediate vicinity of the BH, 
the BLR, traces the gravitational field of the BH. The 
width of the broad emission lines (e.g. H/3) gives the 
velocity scale, while the BLR size is given by the con- 
tinuum luminosity through application of an empirical 
relati on found from reverberation mapping (RM) (e.g., 
Wan del et al.|[l999l: iKaspi et al.ll2000l [2€0a iBeiitz et alj 
(200i). Combining size and velocity gives the BH mass, 

^ Note that electronic tables with all kinematic measurements 
will be presented in the next paper of this series. 
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Fig. 4. — Examples of stellar-velocity dispersion measurements in CaT region, at a distance of ^O-T" from the center. The observed 
spectrum is shown in black, the model in red. The gray-shaded area corresponds to regions that were not included in the fit, due to either 
AGN emission lines or other spurious artifacts. 



assuming a dimensionless coefficient of order unity 
to describe the geometry and kinematics of the BLR 
(sometimes known as the "virial" coefficient). Generally, 
this coefficient is obtained by matching the Mbh-c 
relation of local activ e galactic nuclei (AGNs) to that 
of qui e scent galaxies (lOnken et all 120041 : IGreene fc Hoi 
I2006at iWoo et all |2010( ). Alternatively, the coefhcient 
can be postulated under specific assumptions on the 
geometry and ki nematics of the BL R. We adopt the 
normalizations in iMcGill et al.l (l2008f) . - which are consis- 
tent with those found b y lOnken et al.l ()2004[ ). However, 
since JWoo et al.' (20W) find a slightly different f-factor 
than[Onken et al. (2004), causing a decrease in A/bh by 
0.02 dex, we subtracted 0.02 dex from all BH masses. 

4.3.1. Hp Widths Measurements 

To measure the width of the broad H/3 emission, we 
use the central blue spectrum, extracted with a size of 
1.08" xl". First, underlying broad Fell emission was re- 
moved (if needed) as described in § |3l Then, a stellar 
template was subtracted to correct for stellar absorp- 



tion lines underlying the broad H/3 line in the follow- 
ing manner: We fixed the stellar-velocity dispersion and 
velocity to the values determined in the region ~5050- 
5450A (i.e. the Mgib region that is mostly free of AGN 
emission) and then derived a best fitting-model in the 
region '^4500-5450A, including a polynomial continuum 
and outmasking the H/3 and [OIII] AA 4959,50071 (here- 
after [QUI]) emission lines. The resulting stellar- 
absorption free spectra in the region around H/3 are 
shown in Fig. [12] Finally, we modeled the spectra by a 
combination of (i) a linear continuum, (ii) a Gaussian at 
the location of the narrow H/3 line, (iii) Gauss-Hermite 
polynomials for both [OIII] lines, with a fixed flux ra- 
tio of 1:3 and a fixed wavelength difference, and (iv) 
Gauss-Hermite polynomials for the broad H/3 line. A 
truncated Gauss-Hermite series ()van der Marel fc Franxl 
119931: IMcGill et al.ll200l ) has the advantage (over sym- 
metrical Gaussians) of taking into account asymmetries 
in the line profiles that are often present in the case 
of [OIII] and broad H/3 (Fig. HS]). The coefficients of 
the Hermite polynomials {h^, /i4, etc.) can be derived 
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Fig. 5. — Same as in Fig.lU but for the Mgib region. Note that the spectra end at the rest wavelength 5600A/(l-|-z). (Object 0121-0102 
is not shown here since the S/N is too low to determine (Tspat from this region.) 



by straightforward linear minimization; the center and 
width of the Gaussian are the only two non-linear pa- 
rameters. For [OIII], we allow coefficients up to hi2, for 
II/3 up to /15. 

From the resulting fit, the second moment of the broad 
II/3 ((7H/3) component is measured within a truncated re- 
gion that contains only the broad EE/? line as determined 
interactively for each object (green line in Fig. [T3|. Note 
that the line dispersion is defined as follows. The first 
moment of the line profile is given by 



Ao = y XP{X)dX/ I P{\)dX 
and the second moment is 



X'^P{X)dX/ / P{X)dX 



XI 



The square root is the line dispersion gn a or root-mean- 
squar e (rms) width of the line (see also iPeterson et al.l 
[2OOI . 

We estimate the uncertainty in cth/3 taking into ac- 



count the three main sources of error involved (i) the 
difference between the fit and the data in the region of 
the broad H/3 component (to account for uncertainties 
by asymmetries not fitted by the Gauss-Hermite poly- 
nomials); (ii) the systematic error involved when deter- 
mining the size of the fitting region (which we deter- 
mined empirically to be of the order of 5%), and (iii) 
the statistical error determined by repeated fitting us- 
ing the same fitting parameters (of order 1%). Note 
that due to the very high S/N, the line dispersion in- 
ferred from the Gauss-Hermite polynomial fit is virtu- 
ally indistinguishable from that measured directly from 
the data (on average, fit/data=l±0.02 and at most, the 
difference is 5%). We also compare cth/S with that in- 
ferred from the FHWM assuming a Gaussian profile: 
crH/3/(FHWM/2.355) = 1.11±0.3. The average difference 
of 10%, expected because broad lines are known not to 
be Gaussian in shape, corresponds to a systematic dif- 
ference of 0.04 dex, negligible for an individual object, 
and small compared to the uncertainty on the BH mass 
that we assume (0.4 dex), but potentially a significant 
source of bias for accurate measurements based on large 
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Fig. 6. — Same as in Fig. [4] but in a region redwards of CaHK, excluding the CaHK region due to AGN contamination. (Object 0121-0102 
is not shown here since the S/N is too low to determine a-gpat from this region.) 

to 5100A, assuming a powerlaw of the form oc u" 
with a=-0. 5. (Literature values of a range between - 
0.2 and -1: iMalkan et all 119831: iNeugebauer et all 119871: 
'Cristiani fc Vio' '1990: 'Francis et alJ 119911: iZh eng et alj 
1997; Vandcn Berk et al. 2001; see also [Ben nert ct aQ 
2010 and D. Szathmary et al. 2011, submitted). 

Wc ca lculated BH masse s according to the following 
formula ()McGill et al.ll2008[) : 



samples. 

Fig. [13] shows the fit for all objects. The variety of 
broad H/3 profiles is interesting, with only 6/25 objects 
revealing symmetric line profiles, 8/25 objects having 
more than one peak, and the majority of objects (11/25) 
having asymmetric line profiles, thus showing the need 
for Gauss-Hermite polynomials. While the line profile 
can in principle provide insights into BLR geometry and 
kinematics, this is beyond the scope of the present pa- 
per. Note that the narrow H/3/[0III] A5007A ratio ranges 
between 6-28%, in agreement with o ther studies (e.g., 
iMarziani et all [2001 [Woo et al.|[2006[) . 

4.3.2. Luminosity at SIOOA and BH masses 

We use the SDSS images to simultaneously fit the 
AGN by a point-spread function (PSF) and the host 
galaxy by spheroid and disk (if present). The next 
section and Appendix [X] describe the surface photom- 
etry in detail. The resulting PSF g'-band magnitude is 
corrected for Galactic extinction (subtracting the SDSS 
DR7 "extinction_g" ' column), and then extrapolated 



logA/BH = 7.68 + 2 



log( 



1000 km s- 



0.518 log 



/ XL 



5100 



Vl044 



ergs 



The results are given in Table [4] 

We assume a nominal uncertainty of the BH masses 
measured via the virial method of 0.4 dex. Note that 
we do no t correct for possible effec ts of radiatio n pres- 
sure ( e.g..lMarconi et al.|[200l 120091 see, however, [NetiiB 
[20091 [20m 

4.4. Surface Photometry 

To obtain a host-galaxy free 5100A luminosity (for an 
accurate BH mass measurement) as well as a good es- 
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Fig. 7. — Same as in Fig. [6]for a region including CaHK at ~1.6" from the center. (Object 0121-0102 is not shown here since the S/N is 
too low to determine crgpat from this region.) 



timate of the spheroid effective radius (to measure the 
steUar- velocity dispersion at the effective radius), we per- 
formed surface photometry of archival SDSS DR7 im- 
ages. In our previous papers, we ran the 2D galaxy fit- 
ting program GALFIT (Peng et at 2002^ on high-spatial 
resolution HST images for this purpose. Here, we lack 
space-based images, but - compared to our previous stud- 
ies - the objects are at much lower redshifts (z~0.05 com- 
pared to z~0.4 and 0.6). (Note that the average seeing 
ranges between ~1.2" in the z'-band to ^lA" in the g'- 
band for our sample.) Moreover, the SDSS images come 
in five different filters (u',g',r',i',z'), allowing us to deter- 
mine the host-galaxy properties by simultaneously fitting 
the multiple bands while imposing certain constraints be- 
tween the parameters of each band. Since this is beyond 
the scope of GALFIT, we developed a new image anal- 
ysis code. The advantage of a joint multi-band analysis 
is that it enables us to more easily distinguish between 
the AGN which dominates in the blue bands from the 
host galaxy which is dominant in the redder filters. The 
code is described in detail in Appendix |^ including a 
comparison with GALFIT. The results are summarized 



in Table [S] 

Note that as a sanity check, we checked the clas- 
sic ga, laxy scaling relations (see e.g., iHvde fc Bernardil 
l2009allb[ ). such as fundamental plane (FP) or stellar mass 
vs <Tap,reff and stcUar mass vs rcs^sph- Taking into ac- 
count the small dynamic range and sample size (espe- 
cially when considering elliptical host galaxies only) , the 
results are consistent within the errors. We will show 
and discuss these galaxy scaling relations when present- 
ing the full sample in the upcoming papers of this series. 

4.5. Stellar and Dynamical Spheroid Mass 

Our surface photometry code gives spheroid, disk, and 
total host-galaxy magnitudes for four different SDSS fil- 
ters (g', r', i', z') which can in turn be use d to estimate 
stellar spheroid masses. For this purpose, lAuger et al.l 
(2009) have developed a Bayesian stellar mass estimation 
code that we use here. The code allows informative pri- 
ors to be placed on the age, metallicity, and dust content 
of the galaxy and uses a MCMC sampler to explore the 
full parameter space and quantify degeneracies between 
the stellar population parameters. We use a Chabrier 
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Fig. 8. — Spatially-resolved stellar- velocity dispersions and velocities. In the first and third column, we show spatially-resolved stellar- 
velocity dispersion (upper panel) and velocity curve (lower panel) as derived from CaHK (blue pentagons), Mgib (green triangles), and the 
CaT region (red squares). In the second and fourth column, the corresponding SDSS-DR7 multi-filter image is shown in gray scales. North 
is up, east is to the left; same dimensions as in Fig.[T] with the position of the longslit indicated as black line (corresponding to 20" for 
clarity, but the slit is 175" long). The "-t" and "-" signs indicate the extraction direction of the spectra, corresponding to the x-axis in the 
plots. Note that all figures are on the same gray scale to allow for comparison. 



initial-mass function (IMF). 

Also, with the knowledge of (Tap,roff (as determined 
from the CaT region) and rcff,sph, we can calculate a 
dynamical mass: 

-^^sph,dyn — '^''off .sphCTap.rcff 



with G = gravitationa l constant. For compar ison with 
literature (in particular iMarconi fc Hunt! 12003! ) , we use 
fc = 3. For the same reason, we choose a^n^reS instead of 
Cspat.roff- The results are given in Table [6l 

5. COMPARISON SAMPLES 

For the AfeH-scaling relations, we compile compar- 
ison samples from the literature, including local in- 
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offset from center (srcsec) 
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activ e gala xies i Mar c oni _fc HuntJ [20031 : IHaring fc Ri^ 
120041: f Giiltckin et al.' 2009^ and loca l active galaxies 
(|Greene fc Hq, 2006a : , Woo et al.ll2010t ). Note that while 
for the inactive galaxies, BH masses have been derived 
from direct dynamical measurements, the BH masses for 
active galaxies are calibrated masses either from rever- 
beration mapping or from the virial method. 

5.1. Mbh-ct Relation 

For the Mb h- ct rel ation, we use the data from 
IGiiltekin et al.l (|2009l) (local inactive galaxies). 



Greene fc Hd (l2006al ) (local active galaxies), and 
Woo et all (|2010l ) (local RM AGNs). In all cases, the 
stellar-velocity dispersion measurements correspond to 
luminosity- weighted stellar- velocity dispers i ons w ithin 
a given aperture (Tap. For IGiiltekin et al.l ()2009f ). the 
aperture is typically the effective radius, but as (Tap 
is compiled from the literature, there are also cases 
where it is (Tap,i/8roff or (Tc. However, IGiiltekin et al.l 
lj2009( ) conclude that the systematic differences are small 
compared to other systematic errors. The BH masses 
were determined from direct dynamical measurements 
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Fig. 10. — The same as in Fig. |8] 

(stellar or gaseous kinematics or masers). In total, 
data are a vailable for 49 objects with z < 0.04. For 
I Greene fc~ H o (2006a), (Tap, 1.5" was determined from the 
fiber-based SDSS spectra and is thus within an aperture 
of 1.5" radius. From their s ample of 56 Seyfert-1 galaxies 
with z < 0.1, IWoo et all ([2008^ chose a sub-sample of 
49 objects for which they measured BH mass using the 
line dispersion of H/3 and the Ha luminosit y consistently 
calib rated to our BH mass measurements (|McGill et al.l 
I2008f) . These 49 objects have 5 objects in common with 
our sample, so we use our results instead, leaving us with 
44 local SDSS AGNs. (We will perform a comparison 
for the objects in common to both samples once we have 
our full sample available for which we expect to have a 
total of '^20 objects in common.) Finally, we include 
24 local Seyfert-1 galaxies (z <0.09) for which the BH 
mass has been determined via RM (IWoo et al.l[2010l) . 
For these objects, the a^p measurements were measured 
within an aperture of typically '^l"xl.5" to ~1.5"x3". 

To compare our results with these literature data which 
all consist of luminosity- weighted stellar- velocity disper- 
sions within some aperture, we use Cap.i.s" as determined 
from the CaT region (which is considered the bench- 
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mark). (Note that choosing r^ff instead as aperture size 
does not change the results within the errors; see Ta- 
ble H) 



5.2. Mun-L., 



BH--f^sph 



d A/BH-Afsph Relations 



For the MBH--^ sDh relation, we a gain u se the local inac- 
tive sample from 'Gulte kin et al.l ( 2009f ). here limited to 
35 elliptical and SO galaxies with a reliable spheroid-disk 
decomposition. 

For the Mrh- MsdIi.i^ relation, w e use t he J, H, and K 
magnitudes from lMarconi fc Hunt! ()2003l ) for their group 
1 (i.e. with secure BH masses and reliable spheroid lumi- 
nosities) to calculate stellar masses in the same way we 
calculated our stellar masse s. Also, we updated t he BH 
masses using those listed in iGiiltekin et al.l (|2009f ). This 
leaves us with a sample of 18 objects. 

Finally, for the MeH-Afsph.dyn relation, we com- 
pile local ina c tive galaxies using BH masses from 
IGiiltekin et al.l ()2009D . and calculate dy i iamic al masses 
using ri-eff.sph giv en bv iMarconi fc HuntI (|2003D and CTap 
measurements bv IGiiltekin et al.l (|2009| ). 

6. RESULTS AND DISCUSSION 
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disLaiice from ceriLer (arcsec) 



Fig. 11. — Stellar-velocity dispersion as der ived from the CaHK region (blue pentagons), the Mgib region (green triangles), and the CaT 
region (red squares). In contrast to Fig. 181101 the unresolved data is shown here, i.e. using aperture spectra with aperture widths increasing 
in one pixel steps (corresponding to arcseconds as indicated on the x-axis). 



Here, we describe and discuss our results. After a brief 
section on host-galaxy morphologies, merger rates and 
rotation curves, we focus on the different methods to 
derive stellar-velocity dispersions and perform a quanti- 
tative comparison. Finally, we present the different BH 
mass scaling relations and compare our results to litera- 
ture data. Since the aim of this paper is to outline the 
methodology and present the results of our pilot study, 
we postpone any detailed quantitative conclusions to the 
upcoming papers, once the full sample is available. 

6.1. Host-Galaxy Morphologies, Merger Rates, And 
Rotation Curves 

Using the multi-filter SDSS images shown in Fig. [Tl 
we can determine the overall host-galaxy morphologies. 
Given the low spatial resolutions, we divide the sample 
into three categories: ellipticals (E), SO/a, and spirals 
later than Sa (S). 11/25 objects can then be classified 
as S, 9/25 as E, and 5/25 as SO/a. One object with a 
spiral-like host galaxy morphology is clearly undergoing 
a merger (0206-0017). 1419-1-0754 shows irregular struc- 



ture differing from normal spiral arms and might also be 
in the process of merging. 

The fraction of ellipticals (36±11%) is somewhat 
higher than expected, given that these are (almost all 
radio-quiet) Seyfert galaxies for which the majority has 
typi cally been found to r eside in spirals or SO (~80%, 
e.g. IHunt fc Malkanlll999l and references therein). How- 
ever, due to the low-resolution ground-based imaging 
data, and the small number statistics, there is still some 
uncertainty in this classification, with some objects po- 
tentially falling in the neighboring category. Also, we 
cannot exclude that in a few cases, the images are too 
shallow to reveal the presence of the disk. Indeed, the 
majority of objects (~13/22) show rotation curves with a 
maximum velocity between 100 and 200 kms~^. Also the 
object with a clear merger signature (0206-0017) shows a 
prominent rotation curve with a maximum of 200 km 
rotational velocity, hinting at a spiral galaxy experienc- 
ing a merger event. Both, the variety of host-galaxy mor- 
phologies, in particular with a substantial fraction of host 
galaxies having prominent spiral arms and disk, Hubble 
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Fig. 12. — Spectra around the broad H/3 emission before (upper spectra; red) and after subtraction of stellar absorption and continuum 
(lower spectra; black). 



types Sa and later (44±5%), and the rotation curves un- 
derscore that their kinematic structure is complex and 
indeed spatially-resolved information is necessary. 

The merger rate (0.06±0.02) is lower than for 
our higher redshift objects (0.29±0.1 at z ~ 0.4; 
I Bennert et al.l |2010( ) and more com parable to iiiactive 
galaxies in t he local Univers e (e.g.. IPatton et all 120021 
see, however. iTal et al"1l2OO90 . The merger rate is likely 
to be a function of galaxy mass, with higher (major) 
merg er rates for higher mass objects (e.g.. lHopkins et all 
120101) . Indeed, the local sample studied here has, 
on average, lower host-galaxy luminosities (25 objects; 
logI/host;L0 = 10.33±0.05; rms scatter: 0.29) than the 
one at z ~ 0.4 (34 objects; logLho s t;Lo = 10.54±0.03; 
rms scatter: 0.18; iBennert et al.l ()2010D ). indicating 
lower mass objects. However, we suffer from low num- 
ber statistics and will get back to this discussion, once 
we have analyzed our full sample of '-^lOO local Seyfert-1 
galaxies. 

Note that the image quality does not allow us to de- 
termine the fraction of bars present in the host galaxies 



to st udy the effect of bars on the Mbh-c scaling relation 
(e.g.. iGraham et"ani2010L and references therein). 

6.2. Stellar- Velocity Dispersions 

We here compare the various stellar-velocity disper- 
sions; first spatially-resolved vs aperture stellar-velocity 
dispersions, then the results from three different spectral 
regions. 

6.2.1. Spatially-resolved vs. aperture stellar-veloctty 
dispersions 



Figs. ISlfTOl show the spatially-resolved velocity disper- 
sions for the sample. We are tracing the velocity disper- 
sion for the central 2-6" radius, depending on the object. 
For the majority of objects (~17/22), the overall behav- 
ior of crspat,CaT cau bc described as decreasing from a 
central value of ~130-200 kms^^ to a value of '--^50-100 
kms^^ in the outer parts. For the remaining objects 
(5/22), CTspat.CaT is roughly constant with radius, within 
the errors. 

A different picture emerges when looking at the aper- 
ture stellar-velocity dispersions in Fig. [11] Here, cat- 
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Fig. 13. — Determination of the second moment of the broad H/3 emission. In addition to a continuum, the narrow H/3 component was 
fitted by a Gaussian while the broad H/3 and both [OIII] lines were fitted by Gaussian-Hermite polynomials (see text for details). In some 
cases, strong broad and narrow Hell emission overlaps with the broad H/3 component and was fitted additionally with a broad and narrow 
Gaussian. The observed spectrum is shown in black, the total fit in red and the residual in blue (only in the region that was fitted). The 
region of the fit in which the second moment of the broad H/3 component was determined is shown as green line. Note the variety of H/3 
profiles. For 1535+5754, prominent [Felll] A4658A and [ArlV] AA4711,4740A emission is present in the spectrum and each of these emission 
lines was fitted by a Gaussian component. 



egorizing the overall behavior of iJap.caT as a function 
of "radius" (which is here the increase in width of the 
aperture spectrum), splits the sample into three cate- 
gories, with the majority of objects (9/22) showing a 
constant (Tap,CaT, and the rest distributes with roughly 
even numbers on either decreasing (6/22) or increasing 
(7/22) fTap.CaT with radius. Looking at individual ob- 
jects, 6 objects have a decreasing (Tap,caT with radius in 
the spatially-resolved spectra but shift towards an ap- 
parently increasing crap,caT in the aperture spectra due 
to the rotational support (as reflected in the rotation 
curve). However, overall, the aperture dispersions change 
only slowly with rad ius, a.s has already been no ted by 
ICapellari et"all (|1996D and lGebhardt et al.l (|2OO(]0 for in- 
active galaxies. 

We can make a more direct comparison between 
o'spat,reff determined from the spatially-resolved spectra 
with that determined from the aperture spectra (lap^roff, 



choosing the stellar-velocity dispersion determined from 
the CaT region (see also Table [3]). We divide the sam- 
ple into three sub-categories according to the host-galaxy 
morphology (Elliptical, SO/a, Spiral) and additionally 
distinguish between host galaxies seen face-on and edge- 
on. Fig. [14] shows the result both as function of effective 
radius and (Xspat.reff- The general trend conforms to our 
expectations: If a spiral galaxy is seen edge-on, the ro- 
tation component can bias dap^roff towards higher values 
and thus, "triangles" are expected to lie below the unity 
line. However, since the disk is kinematically cold, it 
can also result in the opposite effect, i.e. biassing crap,roff 
towards smaller values, if the disk is seen face on, so "cir- 
cles" are expected to lie above the unity line. On average, 
face-on spiral host galaxies objects have cTgpat^rcff/o'rGff.ap 
= 1.02±0.05 (rms scatter: 0.24) and edge-on objects have 
o'spat,roff/o'reff,ap = 0.90±0.03 (rms Scatter: 0.16). If we 
calculate the average for the whole sample (all morpholo- 



18 



Bennert et al. 



TABLE 3 

Stellar- Velocity Dispersion Measurements 



Object 


'-'"spat,rcff 


•^ap.rcff 


Cap, 1.5" 


'''cff .sph 




(kms-i) 


(kms-i) 


(kms^-"-) 


(kpc) 


(1) 


(2) 


(3) 


(4) 


(5) 


0121-0102 


104+6 


89+10 


89+8 


1.54 


U20u-0U1 / 


222+2 


200+9 


213+8 


3.07 




159+6 


168+15 


179+11 


1.29 


U5U2-|-v5lU4 


95+2 


128+5 


122+6 


3.17 


Oo4d-|-2o22 


130+3 


205+11 


190+11 


4.48 


1042+0414 


60+1 


85+4 


86+4 


2.76 


1049+2451 


138+4 


189+15 


181+11 


2.78 


llUl + 1102 


113+2 


167+10 


178+9 


4.02 


lllo+412o 


104+2 


110+5 


130+4 


3.40 


11/1,1 1 o/^cro 

1144+^000 


173+4 


154+12 


158+13 


1.26 


1 oi n 1 QQon 
iziU+oozU 




1 1 
lo / itD 




n /I Q 
U.4o 


1250-0249 


112+2 


102+8 


106+9 


2.23 


1323+2701 


113+3 


119+8 


121+10 


1.53 


1405-0259 


119+2 


119+5 


117+5 


1.24 


1419-1-0754 


181+2 


198+8 


194+9 


2.14 


1434-1-4839 


115+2 


118+5 


126+6 


2.16 


1535-1-5754 


99+1 


110+5 


114+5 


2.09 


1545-1-1709 


157+2 


182+6 


182+6 


1.73 


1554-1-3238 


136+3 


152+5 


152+5 


0.66 


1605-1-3305 


134+6 


97+23 


115+12 


0.79 


1606+3324 


130+3 


163+7 


163+7 


1.62 


1611+5211 


95+1 


122+3 


121+2 


2.41 



Note. — Col. (1): Target ID (based on RA and DEC). 
Col. (2): Spatially-resolved stellar- velocity dispersion within 
spheroid effcetive radius from CaT region. Random errors 
are given, while systematic errors are of order 7-15%. Col. 
(3): Aperture stellar- velocity dispersion within spheroid ef- 
fective radius from CaT region. Random errors are given, 
while systematic errors are of order 7-15%. Col. (4): Aper- 
ture stellar- velocity dispersion within 1.5" (to compare with 
SDSS fiber data) from CaT region. Random errors are 
given, while systematic errors are of order 7-15%. Col. (5): 
Spheroid effective radius (in kpc; semi-major axis) from sur- 
face photometry (see Appendix^pand Table[5] fiducial error 
0.04 dex). 

gies and orientations), crspat,refT/cap,reff = 0.93±0.04 (rms 
scatter: 0.2). There is no obvious trend with either bulge 
mass, BH mass, or effective radius. 

To compare the steUar-velocity dispersion measure- 
ments with what would be measured from SDSS fiber 
spectra, we use crap,i.5"- For aspat,refT/a'ap,i.5" , the same 
trend persists as for (Tspat,reff / fap.reff , showing that choos- 
ing 1.5" instead of effective radius does not have a large 
effect on the resulting stellar-velocity dispersion. This 
can be attributed both to the luminosity weighting with 
a stee p central surface-brightness profile ( dc Vaucoulcurs 
119481) and to the fact that the average effective radius of 
our sample is close to 1.5" (2.6±0.07"; rms scatter: 1.8). 

For objects at higher redshift, the effect can be more 
pronounced as different sizes are involved. Consider- 
ing our z~0.4 Seyfert-1 sample for which we study 
the evolution of the MeH-scaling rela.tions (i Trcu ct al. 
2OOIIW00 et al.l[200l[Treu et al.ll2007l: Imoet al, 2008; 
Bennert et al.ll2010[ ). the typical extraction window to 
determine dap is 1 square-arcseconds , with 1" corre- 
sponding to 5 kpc at that redshift ()Woo et al.l 120081 
J.-H. Woo et al. 2011, in preparation). This is a fac- 
tor of ~2.8 larger than the actual effective radius de- 
termine d from surface-brigh tness photometry for these 
objects (|Bennert et al.|[2010l excluding objects with only 
upper limits of the spheroid radius) . If we make another 
comparison, crspat,refT/o-ap,2.8 refr = 0.91±0.05 (rms scat- 
ter: 0.2): On average, aperture spectra can overestimate 



TABLE 4 
BH Mass Measurements 



Object 
(1) 


(kms"-"-) 
(2) 


-^^5100 
(lO*-* ergs-1) 
(3) 


log 


A/bh/A/q 
(4) 


U121-U1U2 


lol ( ±DD 


0.24 




7.58 


OzUb-0017 


1 nn 1 1 1 nn 

1991+iuo 


0.61 




8.15 


rxocTO r\c;oo 

\j6o6-\j'oa6 


1d94±o5 


0.09 




7.58 


U0O2+0IO4 


14/2± /4 


0.05 




7.31 


r\Q> A C 1 ocoo 

Uo4d+2522 


A A '7 I OOT 

4547+227 


0.24 




8.66 


1042+0414 


1 1 ^0 

1252+63 


0.04 




7.12 


10/10 1 liner 

1043+1105 


1 n 1 n 1 n£? 

lyioiyb 


0.16 




7.81 


in/io 1 O/ic^i 


00 c; o_L 1 1 Q 
22o2±llv> 


0.18 




7.98 


1 1 ni 1 1 1 no 
iiUi+liOz 


2yU0±14o 


0.05 




7.91 


1 ,1100 
1116+4123 


1 r\cr 1 1 n cr 

2105+105 


0.03 




7.51 


11/1/1 1 Q^Jt^Q 

1144+3653 


C C 1 _L TOO 

2ool±12o 


0.01 




7.50 


1010 1 ooor\ 

1210+3820 


'7'y 1 1 1 n 

2377+119 


0.03 




7.64 


locro r\oio 

1250-0249 


oo'70 1 1 1 n 

2378+119 


0.06 




7.79 


1323+2701 


2133+107 


0.02 




7.12 


1355+3834 


3110+156 


0.09 




8.10 


1405-0259 


1343+67 


0.05 




7.22 


1419+0754 


1932+97 


0.08 




7.65 


1434+4839 


1572+79 


0.14 




7.62 


1535+5754 


2019+101 


0.26 




7.97 


1545+1709 


1604+80 


0.06 




7.42 


1554+3238 


1988+99 


0.11 




7.77 


1557+0830 


2019+101 


0.08 




7.69 


1605+3305 


1981+99 


0.23 




7.92 


1606+3324 


1737+87 


0.03 




7.36 


1611+5211 


1843+92 


0.04 




7.49 


Note. — 


Col. (1): Target ID (based on 


RA 


and DEC). 



Col. (2): Second moment of broad H/3. Col. (3): Rest-frame 
luminosity at 510QA determined from SDSS g' band surface 
photometry (fiducial error 0.1 dex). Col. (4): Logarithm of 
BH mass (solar units) (uncertainty of 0.4 dex). 

the spheroid-only stellar- velocity dispersion by 0.03±0.02 
dex. This is attributable to a rotational broadening in 
edge-on objects. Note that the fraction of spiral host 
galaxies in the sample at z^O.4 (~14/34) is comparable 
to the fraction in the local sample studied here (~9/25) 
and thus, such a comparison is straightforward. 

However, the bias introduced by measuring stellar- 
velocity dispersions from aperture spectra, even with 
extraction windows much larger than the effective ra- 
dius, cannot explain the observed offset of z^OA Seyfert- 
1 galaxies from t he MBH-Afsph scaling relation seen in 
IWoo et alJ (12001 : For a given BH mass, aap is too low 
for the high-z Seyfert galaxies - the opposite effect to the 
average bias determined here. CTap can be underestimated 
in case of face-on spiral galaxies with a contribution of 
the dynamically cold disk, but this effect is too small 
(on average less than 0.01 dex when considering face-on 
spirals only; see above). We performed the same com- 
parisons also for the other spectral fit ting regions (sinc e 
the region around Mgib was used in iWoo et aD l2008| ) , 
finding similar results. To conclude, aperture effects can 
introduce a small bias in the a measurements but cannot 
explain the offset seen in the Mb h-c relation for higher- 
redshift objects (jWoo et al.l l200g). 

Another more recent study that benefits from our com- 
parison between stellar-velocity dispersions derived from 
aperture spectra to tho se derived from spati ally-resolved 
spectra is the one by iGreene et al.l (|20100 : They re- 
port that a sample of water megamaser residing in spi- 
ral galaxies in the local Universe (~ 0.01 < z < 0.03) 
for which BH masses were derived directly from the dy- 
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TABLE 5 
Results from Surface Photometry 



Object PSF Spheroid 





g 

(mag) 


r 

(mag) 


1 

(mag) 


z 

(mag) 


g 

(mag) 


r 

(mag) 


1 

(mag) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


0191 0109 


17 


01 


17 


51 


17 
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Note. — Col. (1): Target ID (based on RA and DEC). Col. (2-5): Extinction-corrected g', r', i', and z' PSF magnitudes (with an uncertainty of 
0.5 mag). Col. (6-9): Extinction-corrected g', r', i', and z' spheroid magnitudes (with an uncertainty of 0.2 mag). Col. (10-13): Extinction-corrected g', 
r', i', and z' disk magnitudes (with an uncertainty of 0.2 mag). Col. (14): Logarithm of spheroid luminosity in rest-frame V (solar units). Col. (15): 
Spheroid effective radius (in kpc; semi-major a-xis). 



TABLE 6 

Dynamical And Stellar Masses 



Object 
(1) 


OgMsph,vir/Mo 

(2) 


logM h,*/Mo 
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logMdisk,,/M0 
(4) 


logMhost,*/A4 
(5) 


0121-0102 


9.866 


10.12+0.24 


10.60+0.24 


10.75+0 


23 
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11.17+0 
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Note. — Col. (1): Target ID (based on RA and DEC). Col. (2): Dynamical spheroid mass 
calculated from Tcff.sph and CTap.rcff (determined from CaT region; fiducial error 0.1 dex). Col. 
(3): Stellar spheroid mass (using Chabrier as IMF). Col. (4): Stellar disk mass (if present). 
Col. (5): Stellar host mass (only listed if disk present, i.e. if different from (3)). '^: For three 
objects, dynamical masses could not be determined as <Tap,reff could not be reliably measured. 
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namics of the H2O masers fall below the Mbh-o' relation 
defined by inactive elliptical galaxies. As pointed out by 
the authors, given the nature of the host galaxies of these 
megamasers - early-to-mid-type spirals - a bias of the 
stellar velocity dispersion measurements from aperture 
spectra due to the presence of the disk is expected. In 
principle, a rotational component could bias the stellar- 
velocity dispersion measurements to higher values and 
result in the observed offset. However, out of their eight 
objects, only three are significant outliers (their Figs. 8 
and 9), namely NGC 2273, NGC 6323, and NGC 2960. 
Two of these are seen close to face-on (NGC 2273 and 
NGC 2960; their Fig. 6 and 7) in which case the ef- 
fect of the disk component on the stellar-velocity dis- 
persion measured from aperture spectra cannot explain 
the observed offset. Only for NGC 6323, for which the 
disk is seen close to edge-on (their Fig. 6) could the 
observed offset indeed be due to rotational broadening. 
From our Fig. [T3J we estimate that the stellar velocity 
dispersion measured from aperture spectra can overesti- 
mate that from spatially resolved measurements by up 
to ^-^40% in case of rotational broadening by a disk com- 
ponent seen edge on - large enough to move NGC 6323 
close to the local relation defined by ellipticals. To con- 
clude, while for individual objects, the effect of a disk 
on the derived stellar-velocity dispersion from aperture 
spectra can be significant , it cannot explain the offset 
observed bv I Greene et al.l pOlO.) for their full sample of 
megamasers. 
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Fig. 14. — Comparison between different steffar-velocity disper- 
sion measurements, in aff cases determined from the CaT region. 
Ratio between Cgpatreffi i-e. luminosity- weigfited stellar- velocity 
dispersion within effective radius as determined from spatially- 
resolved spectra and a-^p y-^g as determined from spectra extracted 
with an aperture equal to the effective radius, as function of ef- 
fective radius. The upper panel shows all objects for which the 
host galaxy has been classified as ellipticals, the middle panel ob- 
jects with so/a- type host-galaxy morphology, and the lower panel 
objects hosted by spirals. The triangles (circles) correspond to 
galaxies seen edge-on (face-on). The error bars are omitted for 
clarity, measurement errors on the ratios range between 5-10%. 

6.2.2. Stellar-velocity dispersions from different spectral 

regions 

When comparing stellar-velocity dispersion measure- 
ments from the three different spectral regions, the over- 
all picture is that they agree within the errors. A more 



quantitative comparison is shown in Fig. 1151 For the 
aperture data, three extreme outliers were excluded in 
this figure due to contaminating broad AGN emission 
lines and featureless continuum swamping the blue spec- 
tral region (namely 0353-0623, 1049-1-2451, 1535-1-5754; 
see also Fig.fTT]). These outliers are shown in the one-to- 
one plot of (TcaT vs. CTMgib and (TcaHK within the effective 
radii as open symbols (Fig. ITSb ) . (Note that none of the 
objects was excluded for the spatially-resolved data, ex- 
plaining the higher scatter.) 

In general, for the spatially-resolved stellar-velocity 
dispersions (Fig. [TSk . left panels) both ratios 

fspat.CaT/CTspat.Mglb and O'spat.CaT/o'spatCaHK SCattCr 

at the 20-30% level. The average crspat,CaT/o-spat,Mgib 
ratio shows a slight dependence on radius in the sense 
that measuring (Tspat from the Mgib region in the 
center tends to underpredict the "true" (Tspat (here 
assumed to be crspat,CaT) by on average 5-10% while 
it gets overpredicted in the outerparts by on average 
~10-15%. This trend with radius can be attributed 
to the AGN contamination from emission lines and 
featureless continuum that is only present in the inner 
most spectra, and results not only in an increased error 
in the determined velocity dispersion but also in a 
possible bias. The ratio o-spat,CaT/crspat,CaHK, on the 
other hand, shows no clear trend with radius; generally, 
measuring tJspat from the CaHK region overpredicts the 
"true" CTspat by on average 5-10%. Note that there is no 
obvious trend with galaxy morphology. 

For the average ratio from the aperture spectra 
(Fig. [TSb . right panels), the trend is similar to the 
spatially-resolved stellar-velocity dispersions for small 
apertures. At larger radii, the average ratio for 
both o-ap,CaT/o-ap,CaHK and crap.caT/cap.Mgib approaches 
unity. For crap,caT/o'ap,Mgib this is probably due to the 
fact that the dependence on radius seen in the spatially- 
resolved ratio cancels out. 

For the stellar velocity dispersion within the effective 
radius (Fig. 115b ). all ratios are unity within the errors 
(o-spat.CaT/cspat.CaHK = 1.01±0.06; rms scattcr = 0.28; 
CTspat.CaT/cTspat.Mgib = 1.03±0.04; rms scatter = 0.22; 
CTap.CaT/cTap.Mgib = 0.98±0.03; rms scatter = 0.15), ex- 
cept for the stellar-velocity dispersion measured in the 
CaHK region from aperture spectra that tend to over- 
predict CTap.CaT by On avCragC 6% (CTap,CaT/CTap,CaHK = 

0.94±0.03; rms scatter = 0.12). 

To estimate the effect of the potential bias induced 
when using the MgIb region to determine the stellar- 
velocity dispersion for the study of the evolution of 
the MsH-Msph s caling relation (a s done for our z^OA 
Seyfert-1 sample: IWoo et alll2008D . we take into account 
that typically, an aperture much lar ger than the act ual 
effective radius (a factor of ^^2.8 for lWoo et al.l[200^ is 
used for extraction of the spectra (see above). We find no 
bias (CTap^CaT/cap.Mgib = 0.98±0.04; rms scatter = 0.17). 

The general conclusion we can draw from this compar- 
ison is that while the CaT region is the cleanest region 
to determine stellar-velocity dispersions, both the MgIb 
region, appropriately corrected for Fell emission, and the 
CaHK region, although often swamped by the blue AGN 
powerlaw continuum and strong AGN emission lines, can 
also give accurate results within a few percent, given high 
S /N spect ra. This is an important improvement over the 
study bv , Greene fc Ho. (2006h) who use fiber-based SDSS 
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spectra (i.e. aperture spectra with a radius of 1.5") for a 
sample of 40 type-1 AGNs for a similar comparison but 
find a bias of the order of 20-30% (in both Mgib, not 
corrected for Fell emission, and CaHK). 

Furthermore, spatially-resolved spectra are very help- 
ful as they allow to eliminate the AGN contamination 
(powerlaw continuum and strong emission lines) espe- 
cially prominent in the blue spectra (CaHK and MgIb 
region) when extracting spectra outside of the nucleus. 

6.3. Mbh Scaling Relations 

We can now create four different BH mass scaling re- 
lations, namely Mbh-ct, MBH-^sph, MBH-Msph,*, and 
MBH-Msph.dyn and Compare our results with literature 
data (§[5]). The resulting relations are shown in Figs. [16] 
and 1 171 The distribution of residuals with respect to the 
fiducial local relations (Table[7]) are shown as histograms. 

In Fig. [m we plot (Tap,i.5" from aperture spectra 
within 1.5"as stellar- velocity dispersion of our sample, 
for comparison with literature data, for which all mea- 
surement were derived from apertures spectra (with dif- 
ferent sizes, see §[5l we here use 1.5" to be comparable 
to SDSS fiber spectra). Overall, our sample follows the 
same A/bh-c relation as that of the other active local 
galaxies and also that of inactive galaxies. For the lo- 
cal AGNs with stellar velocity dispersion measurements 
from SDSS fiber spectra (green data points in Fig. [16]), 
iGreene fc~ Ho (2006a) already noticed that these objects 
seem to follow a shallower Mbh-c relation with an ap- 
parent offset at the low-mass end in the sense that the 
stellar velocity dispersion is smaller than expected. The 
same trend may also to be visibl e for the RM AGNs 
(blue data points, iWoo et al]|2010[ ) and our local active 
galaxies (red data points). However, for our data points, 
this trend can be attributed to five objects (0121-0102, 
0846+2522, 1250-0249, 1535-^5754, and 1605+3305) that 
might simply be outliers, with strong AGN contamina- 
tion, especially in the aperture spectra. However, at this 
point, the available dynamic range is too small to distin- 
guish between a real offset /change in slope or simply a 
rising scatter. 

The AfBH-^sph relation indicates that our sample of 
active galaxies resides in host galaxies that are overlu- 
minous compared to the inactive galaxies (on average 
by 0.15+0.08 dex; rms scatter: 0.4). One potential 
bias here could be that, due to the shallow images, we 
are missing the disk contribution and thus overestimat- 
ing the bulge luminosity for objects that we classified 
as ellipticals and fitted by a spheroidal component only 
(see also § 16.11) . However, the distribution of residuals 
with respect to the fiducial relation of inactive galax- 
ies shows that especially host galaxies classified as spi- 
rals contribute to this offset (offset 0.2 dex±0.07 dex; 
rms scatter: 0.37), arguing against such a bias. In fact, 
such an offs et might not be too s urprising for two rea- 
sons: (i) the iGiiltekin et ah! (|2009 f) sample only includes 
ellipticals and SO/a in the luminosity plot and (ii) the 
enhanced luminosity might be due to starformation trig- 
gered from a same event that triggered the AGN. That 
AGNs are often hosted by actively star-forming galax- 
ies has been found in various studi es at different red- 
shifts (e g., [p iuffmann et al. 2003; Jahnke et all 120041: 
iHickox et al1i2009i: iMerloni et al.ii2010lV 

However, we cannot exclude that at least some of the 



spiral galaxies have pseudobulges which are character- 
ized by surface-brightness profiles closer to exponential 
profil es (jKormendv fc KennicuttI l2004t iFisher fc Drorvl 
120081 e.g.,). As discussed in Appendix 1X1 using a Sersic 
index of n = 1 instead of n = 4 for the spheroid compo- 
nent can decrease the spheroid luminosity by ~0.2 dex, 
thus accounting for at least some of the offset. We will 
explore this effect further when analysing the full sample. 

For both the MBii-Msph,i, and MBH-A-isph.dyn relations, 
our objects seem to follow the relations determined by 
the inactive galaxies. 

Note that we have a small sample size and also a small 
dynamic range in the parameters covered and all four 
BH mass scaling relations presented here show a large 
scatter. Thus, we refrain from discussing the results any 
further at this point but will get back to the local BH 
mass scaling relations in more detail when we have the 
full sample available. 

7. SUMMARY 

To create a local baseline of the BH mass-scaling re- 
lations for active galaxies, we selected a sample of ^100 
local (0.02 < z < 0.1) Seyfert-1 galaxies from the SDSS 
(DR6) with Mbh > lO^M©. Ah objects were observed 
with Keck/LRIS, providing us with high-quality longslit 
spectra. These data allow us to determine, for the first 
time, spatially-resolved stellar velocity dispersions. Here, 
we present the methodology and first results of a pilot 
study of 25 objects. The full sample will be presented in 
the forthcoming papers of this series. 

From the Keck spectra, we measure both spatially- 
resolved stellar-velocity dispersion and aperture stellar- 
velocity dispersions in three different spectral regions: 
around CaHK, around MgIb (after subtraction of un- 
derlying broad Fell emission), and around CaT. We 
present a detailed comparison between spatially-resolved 
and aperture stellar- velocity dispersions as well as stellar- 
velocity dispersions from different spectral regions. Also, 
we determine the width of the H/3 emission line (after 
subtraction of broad Fell emission and stellar absorp- 
tion). 

On archival SDSS images (g', r', i', z'), we perform 
surface photometry, using a newly developed code that 
allows a joint multi-band analysis. We determine the 
spheroid effective radius, spheroid luminosity, and the 
host-galaxy free 5100A AGN continuum luminosity. 

Combining the results from spectroscopy and imaging 
allows us to estimate BH masses via the empirically cali- 
brated photo-ionization method from the width of the H/3 
emission line and the host-galaxy free 5100A AGN con- 
tinuum luminosity. The spheroid effective radius is used 
to determine the luminosity- weighted stellar- velocity dis- 
persion within ri.cff,sph- The spheroid luminosities in four 
different bands are used to calculate stellar masses. Also, 
our results allow us to estimate dynamical masses. We 
can thus study four different BH mass scaling relations: 

Mbh-O-, MBH-^sph, MBH--Msph,*, and MBH-Msph,dyn- 

The main results for the pilot study can be summarized 
as follows. 

• The host galaxies show a wide variety of morpholo- 
gies with a significant fraction of spiral galaxies and 
prominent rotation curves. This underscores the 
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Fig. 15. — Comparison of stellar-velocity dispersion as measured from CaT, CaHK, and Mgib region, a: Ratio of stellar velocity 
dispersions measured from CaT to MgIb (lower panels) and CaHK (upper panels) as a function of distance from the center. The mean is 
shown at each location including the rms error (shorter - in x - horizontal bar) and the error on the mean (longer - in x - horizontal bar). 
Note that the rms scatter is due to both intrinsic scatter and measurement errors which contribute at the level of 7-15%. The dashed line 
corresponds to a ratio of 1. From spatially-resolved spectroscopy in the left panels (crgpat; i.e. the centroid of the extracted spec tra f rom 
which the ratio was measured is at the distance from center that is indicated on the x-axis with aperture sizes as discussed in § 14.11 and 
with the measurement on both sides of the nucleus averaged); from aperture spectra in the right panels ((Tap; i.e. the spectra from which 
the ratio was measured are always centered on the nucleus but the width of the extracted aperture increases as indicated on the x-axis). b: 
One-to-one comparison for the luminosity-weighted stellar-velocity dispersion within the effective radius: CaT versus MgIb (lower panels) 
and CaT versus CaHK (upper panels) for spatially-resolved spectra in the left panels (cTspat.rcff ) and for apertu re s pectra (crap. reft) in the 
right panels, respectively. The open symbols in the left panels indicate three outliers that were excluded in Fig. I15h .. 




Fig. 16. — Left panel: Mbh-t relation for our sample (red open pentagons), the local RM AGNs fbl ue: |Woo et ar||2010|. with the 
blue dashed line th eir best fit) , an d a local sample of AGNs for which a was measured from SDSS (green; 'Greene & Ho' 2006a') and BH 
masses taken from IWoo et al.l II200 6-). The black datapoints correspond to inactive local galaxies from Giiltekin ct al. (2009) (with the 
black dashed line their best fit; see text for details). Uap corresponds to the luminosity- weighted stellar velocity dispersions within a given 
aperture, depending on the sample (see text for details). The error on the BH mass for both our sample and the local sample of AGNs 
(green data points) is 0.4 dex and shown as a separate point with err or bar in the legend, t o reduce confusion of data points. Right panel: 
Distribution of residuals with respect to the fiducial local relation of Giiltekin et aTl II2009I V Lower panel: literature data; upper panel: our 
sample (black: full sample; blue: spirals, red: ellipticals, green: SO/a) 
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Fig. 17. — The same as in Fig. 1161 but for the other three A/bh scaling relations, namely AfBH-^sph (l^ft panels), MeH-A^sph,* (middle 
panels), and MBH-A'/sph,dyn (right panels). Left panels: For the A/eH-^'sph relation, we use the local inactive sample from lGiiltekin et al.l 
(|2009( ^. here limited to 35 elliptical and SO galaxies with a reliable spheroid-disk decomposition. Middle panels: For the Men-Afaph,* re- 
lation, stellar masses were calculated from the J, H, an d K magnitudes from Ma rconi &: Hun^ ([2003) for their group 1 (see text for details). 
Also, BH masses were updated using those listed in|Gultekin et al. (2009). Right panels: For the AfBH-A/sph,dyn relation, we compile local 
inactive galaxies using BH masses from Giiltekin et ah) l| 20091 ^. and calculate dynamical masses using rjcft.sph given bv lMarconi &: HuntI 
(|2003i ^ and Cap measurements bv iGiiltekin et al.. (|2009f ). 



TABLE 7 

Fits to the Local Scaling Relations 



Relation 


Sample 


a 


/3 


Scatter 


Reference 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 



log(MBH/A-fQ) = a-f /31og((Tap/200kms-i) RM AGNs 8±0.14 3.55±0.60 0.43±0.08 Woo+10" 

inactive galaxies 8.12±0.08 4.24±0.41 0.44±0.06 Gultekin+09 

log(AfBH/Af0) = a-|-/31og(Lsph,v/lO^^L0,y) inactive galaxies 8.95±0.01 1.11±0.18 0.38±0.09 Gultekin+09 

log(A/'BH/M0) = « -|-/31og(Afaph */M0) ' inactive galaxies -3.34±1.91 1.09±0.18 0.38±0.1 here 

log(MBH/M0) = a-|-/31og(Msphdyn/M0) inactive galaxies -0.98±1.31 0.84±0.12 0.54±0.08 here 



Note. — Relations plotted as dashed lines in Figs. [16} [17] and used as fiducial relation when calculating residuals. Col. (1): 
Scaling relation. Col. (2): Sample used for fitting. (Note that we do not fit our local sample as the scatter is too large.) Col. 
(3): Mean and uncertainty on the best fit intercept. Col. (4): Mean and uncertainty on the best fit slope. Col. (5): Mean and 
uncertainty on the best fit intrinsic scatter. Col. (6 ): References for fit. "here" means determined in this paper independently. 
" Assuming the virial coefficient log f = 0.72±0.10 l lWoo et ani2010D . 



• We find a lower merger rate than for our higher 
redshift study, comparable to inactive galaxies in 
the local Universe. 

• Determining stellar-velocity dispersions from aper- 
ture spectra (such as SDSS fiber spectra or unre- 
solved data from distant galaxies) can be biased, 
depending on the size of the extracted region, AGN 
contamination, and the host-galaxy morphology. 
An overestimation of the stellar velocity dispersion 
from aperture spectra is due to broadening from an 
underlying rotation component (if seen edge-on), 



an underestimation can originate from the contri- 
bution of the dynamically cold disk (if seen face 
on). However, co mparing with the h igher-redshift 
Seyfert-1 sample of lWoo et all (|2008[ ). we find that, 
on average, such a bias is small (<0.03 dex) and, 
moreover, in the opposite direction to explain the 
offset seen in the A/bh-c relation. 

• The CaT region is the cleanest region to determine 
stellar-velocity dispersion in AGN hosts. However, 
it gets shifted out of the optical wavelength regime 
to be used beyond redshifts of z c± 0.1. Alter- 
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natively, both the Mgib region, appropriately cor- 
rected for Fell emission, and the CaHK region, al- 
though often swamped by the blue AGN powerlaw 
continuum and strong AGN emission lines, can also 
give accurate results within a few percent, given 
high S/N spectra. Spatially-resolved data are very 
helpful to eliminate the AGN contamination by ex- 
tracting spectra outside of the nucleus. 

• The BH mass scaling relations of our pilot sample 
agree in slope and scatter with those of other lo- 
cal active galaxies as well as inactive galaxies for a 
canonical choice of the normalization of the virial 
coefficient. 
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APPENDIX 
SURFACE PHOTOMETRY 

The imaging results of this paper rely on a new surface photometry code that we have developed (§ 14. 4p . Here we 
describe the code in detail and discuss the SDSS photometric data and how the code was used to fit surface brightness 
models to these data. 

SPASMOID: A New Surface Photometry Code 

Surface Photometry and Structural Modeling of Imaging Data, or SPASM OID, written by MWA, is an image analysis 
code designed to supersede the functionality of GALFIT ()Peng et al.ll2002l ). The code employs a Bayesian framework 
wherein the various model parameters (e.g., the centroid, total magnitude, effective radius, ellipticity, position angle, 
and/or Sersic index) can be tied together by priors (i.e., constraints) between model components or across different 
filters. For example, it is straightforward to impose a prior that the redder bands have smaller effective radii than the 
bluer bands, that the AGN has blue colors related by a power-law of frequency, or that the relative position between 
the AGN and the bulge is the same in all bands and the offset between these components must be small. The code 
also allows for a different PSF model for extended objects and point sources (this is useful if the different components 
have different colors), and can even use linear combinations of PSFs for the AGN to account for PSF- mismatch, which 
may be particularly important for HST imaging. 

The code implements a MCMC sampler to explore degeneracies between the parameters and provide robust error 
estimates. A set of reduced data images is provided by the user along with variance images, image masks, meta data 
(e.g., the photometric zeropoints and the pixel scale in each image), and a starting guess for the parameters (this can 
be substantially different than the best parameters, although a closer guess to the "true" value leads to more efficient 
sampling). A likelihood function is defined assuming Gaussian uncertainties on the pixel values as described by the 
variance images, e.g., 

(image - modcl)2 1 
logP = 2^ 2^ -2 ^^og2nat^^^^. 

images pixels image 

Priors on all relevant model parameters and - for more complicated models - hyperparameters are defined by the user, 
and the code uses the PyMC python module to explore the posterior. We have tested our code on simulated and real 
data and have compared our results with those derived from GALFIT. We find that, for priors that approximate the 
implicit GALFIT priors, we are able to reproduce the GALFIT results (ij 16. 3p . 

A code as flexible as SPASMOID is ideally suited for space-based data; the spatial resolution and depth of the SDSS 
photometry limits the utility of such a code, however, it still has its advantages. We now describe how the SDSS data 
are prepared and subsequently modeled with SPASMOID. 

Preparation of SDSS images 

First, we determine the magnitude zeropoint (zp) following the recipe described on the SDSS DR7 webpage'': 

zp — —1 ■ (a + 6 * airmass — 2.5 ■ log(exptime)) 

with a = zeropoint count rate and h — extinction coefficient taken from the "tsfield" header keywords for a given field 
and filter, and exptime — 53.91 sec. Then, the sky was subtracted using either the sky value in the image header (if 



http:/ /www.sdss.org/DR7/algorithms/fluxcal.html#counts2mag 
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present) or as determined independently directly from the image (note that the 1000-counts bias was also subtracted). 
We created noise images (in counts) according to 



/ data + sky dark_variance 

noise = i / : 1 

y gam gam^ 

with data images in counts and gain and dark_variance taken from the "tsfield" header keywords. Then, for convolution 
with the PSF of the Sloan telescope optics, we use a Gaussian with the parameters given in the "ObjAll" table of 
SDSS for a given object and filter ( "mRrCcPSF jilter", "mElPSF jiher" , "mE2PSF jiltcr"). The average seeing s for 
the 25 objects was Sg = 1.35±0.15, Sr = 1.24±0.15, s,=1.14±0.15, s^=1.16±0.14. 
We only use the four filters g', r', i', and z', as u' is generally too faint. 

Running GALFIT for comparison 

For comparison, we ran GALFIT on all objects, in a very similar fashion as described in detail in I Bennert et al.l 
(|2010l). In short, we first fit the central AGN component with a PSF, thus determining the center of the system 
which was subsequently fixed to all components. The PSF used was created as a circular Gaussian with an FWHM 
corresponding to the seeing, derived from the parameter "mRrCcPSF Jilter" by ^/(mRrCcPS F_filter/2) * 2.355. We 
then fitted a two component model only, consisting of a PSF and a Ide VaucouleursI ()1948D profile. The starting 
paramete rs for the G ALFIT runs were taken from the SDSS DR7 catalog, i.e. PSF magni tude ("psf-mags " ) and 
Ide VaucouleursI ()1948D magnitude ("deVmag") for the spheroid. The minimum radius of the Ide Vaucouleura () 19481 ) 
profile was set to 2 pixels (^0.8"), i.e. the minimum resolvable size given the seeing. 

Fitting with SPASMOID 

Using the GALFIT results (PSF magnitude, location of PSF, spheroid magnitude, spheroid effective radius, ellipticity 
and position angle) as starting parameters, we assume d the fo llowing AGN/host galaxy fitti ng procedure. 

We fit the host galaxy by either a single de Vauc ouleursI |i948) profile or by a de Vau couleursI () 19481 ) plus an 
exponential profile to account for a disk, while the AGN point source is modeled as the Gaussian described for the 
convolution PSF in the previous section. All of the components are concentric but the centroids may vary between 
bands to account for imperfect registration of the images. We also fix the effective radius to be the same in all 
bands (i.e., our photometry is similar to the modelMag photometry of SDSS) as well as the position angle and axis 
ratio, although these are free to vary between the bulge and disk component. Again, the minimum radius of the 
Ide VaucouleursI (| 19481 ) profile was set to 2 pixels. 

The normalizations of the profiles ~ that is, the magnitudes of each component - are determined by generating 
models given the structural parameters (centroid, effective radius, axis ratio, and position angle) and finding the best 
coefficients of a linear fit of these models to the data. In principle the magnitudes could be free parameters of the 
MCMC sampler, but taking advantage of the linear nature of the fit allows us to very quickly find the 'optimal' 
magnitudes for each proposed set of structural parameters. 

Finally, depending on the images, residuals, and the statistics, we decide whether a given object is best fitted by 
three com ponents (PSF, sp heroid, disk) o r two com ponents (PSF + spheroid). (This procedure is similar to the one 
adopted in iTreu et alll2007l and Bennert et al.ll201ft ) 

We subsequently applied correction for Galactic extinction (subtracting the SDSS DR7 "extinction_filter" ' column). 
The resulting AB magnitudes were transformed to rest-frame optical bands by performing synthetic photometry on 
an early-type galaxy template spectrum, a procedure traditionally referred to as k-correction. The template spectrum 
initially has arbitrary units, and these units were adjusted so that the synthetic observed frame magnitudes match 
the magnitudes from our photometry. We then evaluated the V-band magnitudes at the rest-frame of the template; 
luminosities were determined by correcting for the distance modulus given our adopted cosmology. The errors on 
extinction and rest-frame transformation are a few hundredths of a magnitude. We estimate an uncertainty of <0.05 
mag (using the scatter in 20 single stellar population templates with ages ranging from 2 Gyr to 8 Gyr). Table [5] 
summarizes the results. 

Note that our model assume s that the host galaxies can be best fitted by either a single Ide VaucouleursI () 19481 ) profile 
or by a ide VaucouleursI (|1948f ) plus an exponential profile. However, some of the spiral galaxies may not have classical 
bulges, but pseudobulges which are characterized b y surface-brightness profiles closer to exponential profiles (e.g., 
iKormendv fc KennicuttI 120041 : iFisher fc Drorvl I2008D . To test the most extreme systematic uncertainties in derived 
spheroid and PSF magnitude, we re-ran our models for those 14 objects for which we fit a bulge plus disk component 
using a Sersic index of n = 1 instead of n = 4 for the spheroid component. The results are comparable to what has 
already been observed bv lBennert et al.l (j20To) (e.g. their Fig. 10b): Decreasing n from 4 to 1 decreases the spheroid 
luminosity by on average 0.4 mag and increases the nuclear luminosity by on average 0.6 mag. At the same time, the 
disk luminosity increases by on average 0.1 mag. Thus, the extreme systematic effect would move those objects up in 
BH mass by on average ~0.1 dex (small, but not negligible compared to the assumed error of 0.4 dex) and towards 
lower spheroid luminosities by ~0.2 dex. However, the image quality does not allow to determine the best Sersic index, 
especially since the objects are complex in nature due to the presence of the AGN for which a perfectly matching PSF 
fit cannot always be achieved and can result in degeneracies between PSF and spheroid. 
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Comparison between GALFIT and SPASMOID Photometry 

A comparison between the results from our fit and GALFIT is shown in Figs. [TMTUl for the AGN+spheroid de- 
composition. Overall, the agreement is good, showing that our code performs as expected. In detail, the spheroid 
magnitudes agree to within 0.01±0.02 magnitudes (rms scatter: 0.08) with the largest difference at the faint end. For 
the PSF magnitude, the difference can be larger, due to a fainter PSF compared to the spheroid; on average 0.36±0.01 
(rms scatter: 0.33). We use this comparison for a conservative estimate of our error bars, i.e. spheroid and disk mag- 
nitude of 0.2 mag, PSF magnitude of 0.5 mag, and total host magnitude of 0.1 mag. Note that we apply constraints 
between the different filters which is not the case for GALFIT which is responsible for some of the discrepancies, e.g. in 
effective radius (Fig. [T9| . While setting these constraints is less important for an AGN+spheroid fit, it helps to solve 
the degeneracies for a three component fit. Also, the PSF differs slightly between our code and GALFIT as we use a 
circular Gaussian for GALFIT but allow for an elliptical Gaussian in our code. 
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Fig. 18. — Comparison between the parameters obtained using GALFIT (G) and SPASMOID (S). First row: Difference between 
spheroid magnitude determined using GALFIT and SPASMOID as a function of spheroid magnitude (SPASMOID) for the four different 
SDSS filters. Second row: The same as in the first row, but for the PSF magnitude. Third row: The same as in the second row, but as 
a function of difference between spheroid and PSF magnitude. 
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Fig. 19. — Comparison between the parameters obtained using GALFIT (G) and SPASMOID (S). First row: Ratio of spheroid effective 
radius determined by GALFIT vs. SPASMOID as a function of spheroid effective radius determined by SPASMOID for the four different 
SDSS filters. Note that only for GALFIT is the effective radius different between the four different filters, but the same for all bands in 
SPASMOID. Second row: The same as in the first row, but for axis ratio b/a. Third row: The same as in the second row, but for 
position angle. 



